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I . Introduction 

The first six months of the contract were spent in developing the 
mathematical models to be used la the control system design. A major task 
that was completed was the development of a computer program which takes 
aerodynamic and structural data supplied by NASA for the ARW-2 aircraft and 
converts these data into state space models suitable for use in modern control 
synthesis procedures. This program has the ability to generate reduced order 
models by eliminating selected modes. Reduced order models of inboard and 
outboard control surface actuator dynamics and a second order vertical wind 
gust model was developed. In addition, an analysis of the rigid body motion 
of the ARW-2 was conducted, and it was shown that the deletion of the 
aerodynamic lag states in the rigid body modes resulted in more accurate 
values for the eigenvalues associated with the plunge and pitch modes than 
were obtainable if the lag states were retained. 

The remainder of the report consists of a summary of results in each of 
the areas outlined above. The details are given in Working Papers contained 
in the Appendix. 

II. Actuator/ Control Surface Models (Working Papers 1 & d) 

A. Elevator . (Working Paper 1) The elevator transfer function to be 


P«ge 2 


uaad la a simple first-order lag 


XLe 


20 


j + Zo 


The allowable control surface activity levels of Mach 0.86 «nd ISOOO ft. are 
■^1° and -1 2*’ deflection and ^80^/ a for a 12 ft/s vertical gust. The bandwidth 
of the elevator Is much less than the lowest flexural frequency (118 rad/sec), 
and the gain of the elevator at this frequency Is 0.167 (-15.34 Db) ; 
therefore. It appears that the elfivator will not be effective for flutte.* 
suppression. 

B. Inboard Aileron . (Working Paper 1). An eleventh-order model with 
third-order numerator dynamics was given for the Inboard aileron. The first, 
fourth, and sixth flexure modes appear to be the most Important In modeling 
flutter. (See Section on Results.) The frequency of the sixth mode Is 218 
rsd/s. In this range of frequencies, a fourth-order approximation of the 
actual Inboard aileron .'ransfer function exhibits a maximum error of 0.23 Db 
In gain and -9 In phase angle. It Is proposed to use the fourth order 
approximation given as 
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The activity for the inboard aileron ie lislted to 20‘’down and 10 ^ up 
deflection and 130 '*/e deflection rate for at 12 ft/a gust at the flutter 
condition. 

C. Outboard Aileron . (Working Paper 3) The modified transfer function 
for the outboard aileron is given in Working Paper 5. (Nets the discussion on 
the Outboard Aileron given in Working Paper 1 is based on m earlier model of 
the aileron and should be ignored.) The exact transfer function is seventh 
order with second order numerator dynamics. A third-order approximation gives 
the some response characteristics as the exact model up to frequencies of 300 
rad/s. This third order model is 



l.77i-X/0^ 
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The outboard aileron has a maximum deflection of ±13 and deflection rate of 
740 /sec for a 12 ft/s rms gust at the flutter condition. 

III. Wind Gust Model (Working Paper 2) A second-order vertical wind gust 
model given below is to be used 
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where 


V 


veitical gust velocity 
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L “ characteristic length (2500 ft) 

V • the forward velocity 

L * 

- white noise Input with Intensity^ ( 
cr ^ rns gust velocity 

At the flutter condition of Mach 0.86 and 15000 feet, <7* *12 At 

the gust test condition of Mach 0.7 and 15000 feet ^7“ '• 59 ft/sec. The control 
system should renuce bending moments 30 to 40X at all stations at the gust 
test condition. 

IV. State Space Model (Working Paper 3) 

Modem control design techniques require that the system be modeled In 
state space form as 

a 


where 

Z “ state vector 

U ■ control vector 
c 

'jfl > white noise Input 

The equatlon'i of motion for the flexible aircraft are given In the form 
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where 

X “ Che rigid body (plunge end plrch) und elaeClc mode 
deflections 

U “ control surface deflections 
m ve’^tlcsl wind gust velocity 

(See next Section and Working Paper 4 for details of this aodal.) 

Incorporating the cranrfer functions for Che control surfaces and the 
vertical wing gust with Che aircraft aodcl results in a state vector, Z, 
consisting of (1) rigid body deflection and rates, (2) flexural displacements 
and rates, (3) elevator angular deflection, c4) inboard aileron angular 
deflection, deflection rate, and acceleration, (3) outboard aileron angular 
deflection a .d deflection rate, (6) wind gust velocity and an associated 
variable (see Working Paper 2), and (7) lag states associated with the 
unsCeac'y aerodynamics. The control vector, consists of Che commanded 

inputs to Che elevator, Che inboard aileron, and the outboard aileron. The 
white noise input is the forcing term for the wind gust model. 

Since there are 10 flexural degrees of freedom, two rigid body degrees of 
freedom, a lag state for each degree of freedom and each reduced frequency 
included in Che unsteady aerodynamic model, a first order elevator mor'el, 
fourth and third order aileron models, and a second order gust model, the 
dimensionality of the state vector is 46 even if only one lag state is 
assumed. Thus it is necessary to computerize the manipulations required to 
construct the state space model. Such a program has been written and is 
running on the University of Minnesota Computer System. This program allows 
rigid body and flexure modes to be delected in order to generate a lower order 
model if required. Results obtained from use of this model are described 


later. 
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Unateady Aarodynulc Mcdel (Working Paper 4) 
The flexible aircraft la modeled aa 
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where 


M > generalized nar>8 matrix 
C •• generalized atructural damping matrix 
R “ generalized stiffness matrix 
q “ dynamic pressure 

Q “ matrix of aerodynamic coefficients 
X, U are defined In Section 4. 

The matrix of aerodynamic Influence coefficients, Q (s-J ^ k) , was provided 

c 

by NASA for a range of reduced frequencies. The aerodynamic Influence 
coefficient matrix was approximated by 


+ z -- 

ryy^t ^ 






The error matrix E(s) Is defined as 
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Determination of matrlcea A^, A^, A^, and D^whlch give the beat leaat aquarea 
fit to the data ia a relatively atraightforward problem (aee Working Paper 4) 

The first column of the matrix A^muat be set equal to the first column 
of QfO^Cwhich is zero) in order to reflect the fact that aerodynamic forces due 
to the plunge displacement are zero. The remaining matrices can be determined 
to give the best least aquarea fit to the data once the k^'s are specified. 
The selection of k 'a which result in a "beat" fit is not so straightforward. 
(See Dowell, E. H. "A Simple Method for Converting Frequency Domain 
Aerodynamics to the Time Domain" NASA Tech. Memorandum 81844, Oct 1980.) An 
approach to the selection of the k 's was developed as part of this s :udy and 
appears to yield good results with minimal computational effort. This 
approach depends upon the fact Chat the spectral norm of a matrix equals its 
maximum singular value. Thus Che maximum singular value of E(s), defined as 
0'C^)equals E(s) and is a measure of the size of the error in the approxlmat on 
of Q(s). The procedure is Co (1) arbitrarily select values for the (2) 

let the first colwuri of A ^ equal the first column of Q(G) (which is zero), (3) 
determine Che remaining values of A^, A^, A^, and , a-1, ...n which give 
Che best least squares fit to the data, ( k) calculate Che maximum singular 
value of Che error matrix E(s), (5) vary Che k^'s and repeat Che process until 
this singular value achieves a minimum. Since only a few values of k^ are 
commonly used, a relatively simple search procedure can be used to determine 
Che optimum values for Che 
VI. Results 


A. Rigid Body Modes (Working Paper 4) 

The procedure described above was used to generate a mathematical model 
of Che ARW-2 at the flutter condition of Mach 0.86 and 13000 feet. A single 
aerodynamic lag state was used. The value of k^was varied until Che maximum 
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singular value of tha error matrix was alnlmlxad for a rsducad fraquancy of 
zero. Tha minimization was accompllshad at aaro fraquancy bacausa It was fait 
that the approximation should ba bast at low fraquanclas slues rigid body 
modas wars to ba studlad. Tha rasultlng valua of k^was 0.13 which Is vary 
closa to tha raducad fluttak. fraquancy of O.IS. Tha alamants of the 0 matrix 
suppllad by NASA wars plotted In polar form and ware compared with polar plots 
of the elamants of tha matrix resulting from tha approximate modal. As can 
ba seen In Working Paper 4, tnesa plots are almost Identical, indicating that 
tha approximation Is axtramaly good. 

All flexure modes ware neglected, and the eigenvalues associated with tha 
plunge and pitch modas ware calculatao as follows: 
plunge 6.1146, 1.828 x 10~^ 
pitch -5.1767 + j 7.5543 

If the lag terms were neglected when the eigenvalues were calculated, the 
plunge and pitch eigenvalues were 

-Z -7 

plunge 3.391 X 10 , 7.588 x lO 

pitch -1.1188 ± J 3.5572 

These are not too different from the values given by Boeing of 
plunge -0.0092 ± J .0437 
pitch -1.4274 ± J 2.422 

It appears that Inclusion of the aerodynamic lag terms degrades the accuracy 
of the calculation of the rigid body eigenvalues and that ellmluatlcn of the 
lag states associated with the rigid body modes Improves the accuracy of the 
model. It is difficult at present to explain why this Is so. The aerodynamic 
lag terms represent unsteady aerodynamic effect and since use of quasi-steady 
aerodynamics usually allows the accurate prediction of rigid body eigenvalues. 
It Is not surprising that the eigenvalues calculated by neglecting the lag 
terms are near those given by Boeing. What Is surprising Is that the 
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inclualon of th# lag Canaa haa such an efface on th* alganvaluaa. Th* Ciaa 
cooaCanta aaaociaCad with Cha unaCeadv aarodynaalca ar* laaa chan .01 
while Che dm conaCanC aaaoclaccd wU)t t^ho shore period node is of Che order 
of 0.5 sec and chac of Che plunge aode Is approximeely 125 sec. 

B. Flexure Hodes 

An exaainacion of Che 10 flexure laod'a indicaCed ChaC nodes 2 and 5 ware 
primarily fuselage bending modes, and mode 7 was exclusively a Call mode. 
Therefore, Chese chree modes were noC considered furcher in Che analysis . 
Mode 1 was the firsc wing bending mode, mode 2 was Che second wing bending 
mode, and mode 6 was Che firsC wing lorsior node. These nodes were obviously 
imporcanc and were reCained. Modes 3 and 8 included wing tip bending and 'c 
was felc Chat Ches . modos should also be analyzed furcher. Mode 9 w^a 
primarily wing bending and mode 10 was primarily wing Corsion, and Chese modes 
were also reCained. Thus seven flexure modes, Che 1st, 3rd, 4ch, 6ch, 8ch, 
9ch, and lOCh, were used in Che fluccer malysis. A single lag Cerm wich 
reduced frequency of 0.13 was used. The loci of Che eigenvalues of chose 
modes as velocity <s varied are shown in Fig 1. The data from which Fig. I is 
constructed is given in Table 1. It can be seen Chat Che 1st mode flutters at 
a speed of approximately 9500 in/s at a frequency of approximately 120 rad/s. 
The results are based or. aerodynamic data for a Mach number of 0.86; 
therefore, Che values of Che eigenvalues at low velocities are suspect, 
however, 9500 in/s corresponds Co a Mach number of 0.75 which is reasonably 
close Co 0.86. Fig 1 presents results which are very similar Co Chose given 
by NASA; therefore, it is felt Chat Che mathematical modeling has been done 
correctly. It should be noted chac modes 3 and 8 are very insensitive to 
velocity, indicating that they are primarily vibrational modes not affected by 
aerodynamics. Modes 9 and 10 do however very considerably wich velocity. 

Eigenvalues were calculated using models in which various modes were 
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dalatad. (Saa Ta'ola 2.) It can ba aaan chat avan chough tha laC aoda 
fluttara, a chraa soda aoda) which contalna tha lac, 4ch, and 6ch aodaa la 
raqulrad Co accuracaiy pradlcc flutCar. Thla la not aurprlalng alnca 
claaalcal fluCCar raqulraa both a brndlng and a Coralon aoda. Both Cha lac 
and 4ch aodaa ara banding and cha 6ch la Coralon. It la InCaraaClng Co noCa 
Chat at Mach 0.6, Cha 4ch aoda fluttara. DalaCion of Cha 3rd and 8ch aodaa 
haa alaoac no affacC on cha alganvaluaa of Cha othar aodaa, and Cha 9ch and 
10th aodaa hava vary IlCCla affacC on Cha lowar aodaa. 

VII. Conclualona and Putura Plana 

It la fait Chat wa now hava a valid aodal of Cha aircraft Co ba uaad In 
our fluCCar control aCudlaa. Since Cha elevator haa auch a low bandwidth 
coaparad with cha flutter frequencies, flutter control will ba accoapllshad 
using Che Inboard and outboard ailerons. A aodal consisting of the Isc, 4ch, 
and 6ch aodaa will be used for control sysCea design and a nodal containing 
Che 1st, 4ch, 6ch, 9Ci , and lOch modes will ba used for evaluation. A single 
aerodynamic lag state Is proposed Co be used In both tha design and evaluation 
models. Currently a program Co perform elganspaca design Is being written an 
Is a prograa to Interface the aircraft model with an exlatlng evaluation 
program for determining rms responses Co stochastic wind gusts. 
























Table 1. Variation of Eigenvalues with Velocity at ISOOO Feet 


Dyn- Pressure (Psl) Velocity (In/s) Mode No. Eigenvalue 

0.0036 1000 iO -4.6 ± j 499.3 

9 -2.8 ± J 421.5 

8 -2.3 ± J 395.8 

6 -4.1 ± j 267.3 

4 -2.3 ± J 191.3 

3 -0.7 ± J 136.7 

1 -1.6 ± J 50.7 

0.144 2000 iO -6.7 ± j 498.2 

9 -3.6 ± j 421.6 

8 -2.5 ± J 395.8 

6 -6.8 ± j 264.9 

4 -3.7 ± J 191.2 

3 -0.7 ± j i36.7 

1 -3.0 ± j 51.6 

0.325 3000 10 -8.9 ± j 496.4 

-4.3 ±J 421.7 

8 -2.8 ± J 39J.9 

6 -9.5 ± j 260.9 

4 -5.1 ± j 191.1 

3 -0.7 ± j 136.7 

1 -4.5 ± j 53.3 

0.577 4000 10 -11.0 ±j 493.8 

9 - 5.1 ± J 42: ,9 

8 - 3.1 ± j 396.0 

6 -12.0 ± j 255.4 

4 - 6.6 i: j 190.3 

3 - 0.8 ± J 136.7 

1 - 6.0 ± j 55.9 

0.902 5000 10 -13.1 ± j 490.5 

9 - 5.9 ± j 422.2 

8 - 3.3 ± j 396.2 

6 -14.3 ± j 248.2 

4 - 8.4 ± j 190.0 

3 - 0.8 ± j 136.7 

1 - 7.7 ± J 59.4 

1C -15.1 ± j 486.3 

9 - 6.7 ± 3 422.5 

8 - 3,5 ± j 396.4 

6 -16.12± j 239.7 

4 -10.7 ± j 188.0 

3 - 0.8 ± j 136.8 

1 - 9.5 ± J 64.4 


1.3 


6000 


Table 

1. Variation of 

Eigenvalues with Velocity 

at 15000 Feet 

i.77 

7000 

10 

-17.16± j 

481.3 



9 

- 7.5 ± J 

422.9 



8 

- 3.7 ± J 

396.5 



6 

-16., 9 ± j 

230.5 



4 

-13.9 ± J 

183.2 



3 

- 0.8 ± j 

136.8 



1 

-11.7 ± j 

71.6 

2.31 

8000 

10 

-19.1 ± J 

475.5 



9 

- 8.4 i J 

423.3 



8 

- 3.9 ± j 

396.8 



6 

-16.4 ± J 

223.0 



4 

-18.3 ± J 

171.5 



3 

- 0.9 ± J 

136.9 



1 

-14.2 i: J 

82.9 

2.92 

9000 

10 

-21.2 ± J 

468.9 



9 

' 9.2 ± j 

423.8 



8 

- 4.0 ± J 

397.0 



6 

-15.7 ± J 

218.8 



4 

-24.0 ± j 

140.1 



3 

- 0.6 ± J 

137.4 



1 

-16.7 ± J 

108.6 

3.61 

10000 

10 

-23.0 ± J 

461.2 



9 

-10.1 ± j 

424.4 



8 

- 4.1 i: J 

397.2 



6 

-16.0 ± i 

217.6 



4 

-68.9 ± J 

117.9 



3 

- 0.5 ± J 

136.6 



1 

-*•21.2 ± J 

121.9 

4.37 

11000 

10 

-24.8 ± J 

452.5 



9 

-11.0 ± J 

425.0 



8 

- 4.1 ± J 

397.4 



6 

-16.9 ± j 

218.3 



4 

-97.1 ± J 

108.4 



3 

- 0.7 ± J 

l.'*6.6 



1 

-*•41.6 ± J 

118.1 



Table 2 


Eigenvalues for Combinations of Flexure Modes 


Mode 


1 

4 

1 

6 

I 

4 

3 
I 

6 

4 
1 

6 

4 

3 
1 

8 

6 

4 

3 
1 

10 

9 

6 

4 
1 

9 

8 

6 

4 

3 
1 

10 

9 

8 

6 

4 
3 
1 


No. 

1 Mode Model 

2 Mode Model 

2 Mode Model 

3 Mode Model 

3 Mode Model 

4 Mode Model 

3 Mode Model 

5 Mode Model 

6 Mode Model 


Eigenvalue 


-13.6 

± 

J 

78.6 

-14.7 

± 

J 

180.2 

-15.0 

± 

J 

86.6 

■•■14.9 

± 

J 

134.2 

-65.0 

£ 

j 

121.7 

-14.7 

± 

J 

180.2 

- 0.9 

it 

J 

136.8 

-15.0 

± 

J 

36.5 

-15.7 

± 

J 

216.1 

-94.3 

± 

J 

106.4 

+39.9 

± 

J 

118.0 

-15.7 

± 

J 

216.1 

-94.2 

± 

J 

106.6 

- 0.7 

± 

J 

136.6 

+39.7 

± 

J 

118.0 

- 4.9 

± 

J 

397.2 

-15.7 

± 

J 

216.1 

-94.1 

± 

j 

106.7 

- 0.7 

± 

j 

136.6 

+39.6 

± 

J 

118.7 

-24.3 

± 

J 

453.6 

-10.2 

± 

J 

425.0 

-16.6 

± 

J 

218.3 

-94.2 

± 

J 

108.7 

+39.8 

± 

4 

J 

118.2 

-10.9 

± 

J 

424.9 

- 4.2 

± 

J 

397.4 

-15.7 

± 

J 

216.1 

-93.7 

± 

J 

107.2 

- 0.7 

± 

J 

136.6 

+39.3 

± 

J 

118.2 

-24.8 

± 

J 

452.5 

-11.0 

± 

J 

425.0 

- 4.1 

± 

J 

397.4 

-16.9 

± 

J 

218.3 

-97,1 

± 

J 

108.4 

- 0.7 

± 

J 

136.6 

+41,6 

± 

J 

118.1 
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Working Paper No. 1 
Actuator Models for 
Flutter Control Study 


W.L. Garrard 

Dept, of Aerospace Engineering and Mechanics 
University of Minnesota 
Minneapolis, MN 


Introduction 

The DAST-2 vehicle has three independent control surfaces 
available for flutter and pitch control. These are an elevator, 
an inborad actuator, and an outboard actuator. This paper gives 
the transfer functions for these actuators and suggests some 
lower order model approximations for the ailerons. 

Elevator 

The elevator transfer function is a simple first order lag 


( 1 ) 

Bode plots for this transfer function are given in Figs. 1 and 
2. The allowable rms control surface activity levels are +7“ 
and -12° deflection and +80°/sec for a 12 ft/s vertical gust at 
the flutter condition of Mach 0.86 and 15000 ft. altitude. 


Inboard Aileron 
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The transfer function for the inboard aileron is 

« 7. 2 3^ )( s-n, !o^( ^ ZgoD<i B. iZd 45 5 5) 

^ r K • 




{ 47/5 4- 2.2 i X 3ZZ 6 4 7. 7o7S \ 

^hS.is ^ 8.ZZ1 X f ZS'^d.'^6 + 


ra ^ -hi. I 71 Xn') 

or in factored form i ^ ^ ^ 


1 



'7. Z3 3 X/<3 s g ^ Moo i. I /(j 79 . a ^ y 

S 2 cra\/' 'J 

A * ® ✓ Cs l ^ \ A m tr r^\y ^ . 


ORIGINAL PAGE IS 
OF POOR QUALITY 


3j^ 7 :f • 54*,/. <s) 




Bode plots for this transfer function are shown in Figs. 3 and 
2 6 

4. The terms s + 2800s + 3.126 x 10 in the numerator and 
s^ + 2848.9 S+ 3.183 x 10^ can effectively be canceled. The 
two high frequency factors in the denominator (s + 1424.5 + 
j 1074.2) and (s + 332.7 + j 8461.6) have little effect on the 
dynamic response over the frequency range of interest in the 
control problem of 1 to 500 rad/s. Also the s + 610 term in 
the numerator comes airly close to canceling the s + 455.5 
term in the denominator and if there are no numerator dynamics 
the state space representation is simplified somewhat. All 
of these assumptions result in a sixth order model 



1.3705 X /O 
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The Bode plots over the frequency range from 10 to 1000 rad/s 
are given in Figs. 5 and 6. It can be seen that the frequency 
response for the exact transfer function and the sixth order 
approximation are very close over this range of frequencies. 

The comparison is presented in detail in Table 1. The flutter 
frequency is about 150 rad/s and the correspondence between 
the sixth order and exact transfer function is very good near 
this frequency. At very high frequencies the exact transfer 
function has a phase shift of -180° and a slope of 40 DB/Decade 
greater than the sixth order approximation. 

A fourth order approximation can be generated by elimi- 
nating the dynamics associated with the factor s + 144 + j895.5. 
This is somewhat difficult to justify as we retain the factor 
s + 161 + j825.8 which is about the same frequency. The result- 
ing fourth order transfer function is 

339. SY 34 SZ5.S ') 
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3ode plots for this transfer function are given in Figs. 7 
and 8 The fourth-order model gives a reasonably good approxi- 
mation of the exact transfer function up to about 300 rad/sec 
but deteriorates rapidly at higher frequencies. The numerical 
details of the comparison are given in Table 1. 

The final approximation considered is second order and is 
given as 

II ' / — 8o z'' 

33£.S 53^.a') 


As can be seen from Figs. 9 and 10 and from Table 1, the accur- 
acy of this model deteriorates fairly rapidly for frequencies 
above about 175 rad/s. 


The activity for the inboard aileron is limited to -10“ 
and +20“ deflection and to/30“/sec deflection rate for a 12 ft/s 
vertical gust at the flutter condition. 


Outboard Aileron 


zU. 


The transfer function for the outboard aileron is 
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Bode plots for this transfer function are given in Figs. 11 and 
12 for frequencies from 0.1 to 1000 rad/s and in Figs. 13 and 14 
for frequencies from 10 to 1000 rad/s. If we consider only fre- 
quencies of the same order as the structural frequencies, a 
fourth order model results. This is given by 


U 




Jl 




I.5S/ 


X /o 


n 


+ szi.:^i'^3Zs.2 Yei Z3s. 


C ± \ 554 '^ 
0 


3 



Bode plots for this transfer function are given in Figs. 

15 and 16. It can be seen that this gives a very good approxi- 
mation of the exact transf ?r function over the frequency range 
of interest although at very high frequencies the exact model 
will exhibit a phase shift of -360® and a slope of -80 DB/Decade 
greater than the fourth order model. 


A second order model given as 


JU 



J 7JS^ 


was also examined and Bode plots for this model are shown in 
Figs. 17 and 18. The approximation resulting from use of this 
model is considerably worse than given by the fourth order model 

The outboard aileron has a rms deflection limit of +15® 
and a deflection rate limit of 710®/sec for a 12 ft/sec vertical 
gust at the flutter flight condition. 

Conclusions and Recommendations 

Initial controller designs will incorporate rigid body and 
first, fourth, and sixth elastic modes. Since the sixth mode 
has a frequency of 225 rad/s at the flutter condition, it is 
felt that the fourth order approximation of the inboard actua- 
tor whould be adequate even though the sixth order gives a 
considerably better approximation of the actuator dynamics for 
frequencies greater than 350 rad/sec. The fourth order model 
of the outboard actuator gives a very good approximation of 
the exact transfer function over the frequency range of interest 
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Table 1: Comparison of Exact Transfer Function for 

Inboard Aileron and Various Lower Order Approximations 


— 

Freq ( rad/s ) 

Gain (Db) 

Phase (deg) 

A Gain (Db) 

A Phase (Deg) 

Exact 

100 

0 

-27 

— 

— 


178 

0.5 

-45 

— 

— 

1 

316 

0.75 

-81 

-- 

-- 


562 

2 . 0 

-162 

— 

-- 

' 

1000 

0 

-396 

— 

-- 

1 

830 

6 . 6 (peak) 

— 


— 

2nd Order 

100 

0 

-18 

0 

-9 

Approx . 

178 

0 

-27 

.5 

-18 


316 

-.075 

-54 

1.5 

-27 

1 

562 

-4.5 

-108 

6.5 

-54 

1 

1000 

-13.0 

-135 

13.0 

-261 

4th Order 

100 

0 

-18 

0 

-9 

Approx . 

178 

0.25 

-36 

0.25 

-9 


316 

0 . 5 

-62 

0.25 

-19 

1 

562 

-0.5 

-135 

2.5 

-27 

1 

1000 

-7.5 

-270 

7.5 

-126 

1 

1 

298 

0 . 5 (peak) 

— 

6.1 

— 

1 6th Order 

100 

0 . 25 

-23 

-.25 

-6 

I Approx, 

178 

0.5 

-40 

0 

-5 


316 

1.25 

-72 

-.5 

-9 


562 

3.0 

-145 

-2.0 

-23 

i 

1000 

4.5 

-300 

-4.5 

-36 

1 

830 

8 . 3 (peak) 

— 

-1.7 
















Table 2: Comparison of Exact Transfer Function for 

Outboard Aileron and Various Lower Order Approximations 
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Working Paper No. 2 

Gust Model for 
Flutter Control Study 

W.L. Garrard 

Dept, of Aerospace Engineering and Mechanics 
University of Minnesota 
Minneapolis, MN 


Introduction 

A second order vertical wind gust model is described. The 
vertical gust is to have a 12 ft/sec rms value at the flutter 
condition of Mach 0.86 and 15000 feet and a 59 ft/sec rms value, 
for the gust test condition of Mach 0.7 and 15000 ft. A 30 to 40 
percent reduction in bending moment at all stations is desired 
at this condition compared with the uncontrolled aircraft. 

Gust Transfer Function 

The transfer function for a second order gust model is given 
as 

f3~ £ ) 

(( I ^ ^ 

where 6 = the vertical gust velocity 

£ = the characteristic length (2500 ft in this case) 

V = the forward velocity 
n = a white noise input 

If we assume s is the spectral of the white noise the expected 
2 

value of 5 is given as 




2 


If constant, this integral can be evaluated from standard 

tables [1] as 


Sts'] 



2 2 

Now we want E[6 ] to equal a , the specified rms value of the gust; 
therefore , 


Now 







J- [ 

2TT J 


od 




where R^(t) is the autocorrelation function of n and since n is 
white noise 


- -A- SCr ) 

ii ^ ^ 


where 6 is the Dirac delta function. Then 



therefore the intensity of the white noise is 



At flutter, Mach = 0.88, h = 15000 ft, v = 908.8 ft/s, Z/v 

/Z 75 a *■ > J 

and 




2.75 


A = 396.00 


At the gust test condition, v = 739.7 ft/s, Z/v = 3.38 
4/0 i 3.3S 5 -^ » / 


and 


A 


11765.78 


3 


terns . 


Both trinsfer functions represent critically damped sys- 
The natural frequency at the flutter condition is 


co^ = 0.363 rad/sec (flutter condition) 
and at the gust test condition is 

Cot^~ 0.296 rad/sec (gust test condition) 

Bode plots for the flutter condition are shown in Figs. 1 and 
2 and for the gust test condition in Figs. 3 and 4, 


State Space Representation 


Since the gust models contain numerator dynamics a little 
extra wor)c is required to put them in state variable form. We 
can accomplish this by using the bloc)c diagram shown in Fig. 5. 
The transfer function is 



By inspection 

h,-- 7U3 -V 

and 

n, : - 

Now from the bloc)c diagram we have the equations of motion in 
state space form as 

g * 0.28S 

For M = 0.86 and 15000 feet 

J X f *■ 0./04X 

r ~ 0.7Z7 yt - O. J3Z S ~ O. OS(^3 % 
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and for M = 0.7 and 15000 feet 


h ~ ^ 0-0^ 

t - O.OS7s^ - O. CX3SS% 


Reference 

1, Crandall, S.H. and Mark, W.D., "Random Vibration" 
Press, 1963, page 72. 
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Working Paper No. 3 

State Space Models for Flutter Control Study 



Working Paper No. 3 

State Space Models for Flutter Control Study 
B. S. Liebst 

Depart-inent of Aerospace Engineering and Mechanics 
University of Minnesota 


Introduction 

The structural, aerodynamic, and actuator models proposed 
for use in the analysis of the DAST-2 vehicle have earlier 
been presented in the frequency domain (see working papers No. 1 
and No. 2) . This paper gives the corresponding state space 
representation of these models for use in the modal control system 
design . 

Second Order Stale Equation 

The following is the proposed second order form of the 
state equations of motion relating vertical gusts and control 
surface deflections to the structural response of the vehicle. 

X = modal coordinates and rigid body modes 

u = control surface deflections 

6 = vertical gust velocity 

V = forward velocity 

q = 

c = mean aerodynamic chord 
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Actuator Model 

The following is proposed frequency and corresponding state 
space models for the three actuators of the vehicle. See working 
paper No. 1 for details of the transfer function representations. 


I . Elevato r 

Ut. o i 
Ut '■0 

therefore , 

tTc. + 2.0 Me. 


2.0 

5 + 2. o 
2-0 Uc. 


II . Inboard Aileron 

II 

LA ^ ^ /. (^ / 4 XI 0 

UiU) c 5 2: j 8)Cf yw / t j 8 ;) 
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.... ^ , 
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III. Outboard Ailerot* 
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Gust Model 


The following is the white noise model for the vertical gust 
state. See working paper No. 2 for details of this derivation. 

6 vertical gust velocity 
2 = intermediate gust state 

n = zero mean white gaussian noise 

i = characteristic gust length 
2 

a = rms value of the gust 



I -- 
« 

S - 

where 


First Order State Equation 

Now that the structural, actuator, and gust models have 
been developed in the state space, they can all be adjoined to 
form one first order state equation. Neglecting the S6 and Pu 
terms of the second order equations of motion the resulting first 
order state equation is 


X = Ax + Bu + w 

where 



Ur- - 
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in 
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NOMENCLATURE : 

M: Generalized mass matrix 

k: Generalized stiffness matrix 

C; Generalized damping matrix 

x; Vector of generalized coordinates for rigid and elastic modes and 

control surface deflections 

Vector of generalized coordinates for gust Inputs 
Q(s) : Unsteady aerodynamic Influence coefficient matrix 

li); frequency In radians per second 

c: reference chord length 

v; velocity of the vehicle 

q: dynamic pressure 

^ : density of air 

II: Mach number 

H: Altitude 

: laplace operator 

a [a] : maximum singular value of matrix A 
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INTRODUCTION 


For the design of active control systems for the suppression of aero- 
dynamic flutter. It Is necessary to first obtain the state-space represen- 
tation of the equations of motion for the flexible vehicle. This is so 
because at present all modern control design techniques are based on the 
availability of a state-space modal. 

In the study of flexible vehicles, the unsteady aerodynamic forces 
and moments are evaluated at various reduced frequencies by the use of 
some type of finite element computational procedure. Our objective is to 
obtain a model of tho unsteady aerodynamic forces and moments in a form 
which can be Incorporated into the structural equations of motion for the 
aircraft so as to get a suitable state-space representation of the vehicle 
dynamics. 

This paper discusses the procedure for obtaining such a model. The 
following discussion presents an approximation of the unsteady aerodynamics 
by a rational polynomial and is based on the approximation first suggested 
by R.T. Jones [ll. The details of obtaining such an approximation using a 
least squares fit over the range of frequencies available is then presen- 
ted . 

A model for the unsteady aerodynamics of the DAST ARW-2 aircraft, 
the data for which was provided by NASA, was developed using the proce- 
dure discussed herein. Using this model, the rigid body motion of the 
aircraft was studied and the results compared with those obtained by 
NASA. 

The programs developed to obtain the least squares approximation 
and do the rigid body analysis are appended for thoroughness. 
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MODEL FOR UNSTEADY AERODYNAMICS: 

The general equations of motion for a flexible vehicle at a given 
Mach number are [2]: 

.a ^ /• y 4- k]^ ^ ^ 


0 ) 




If a simple harmonic motion Is assumed, the matrix of aerodynamic 
influence coefficients Q(/l - k) , where k is the reduced frequency 
given by k is calculated using finite difference procedures at a 


finite number of reduced frequencies k^^, i » 1, 2,..n. 
matrix can be approximated by a polynomial in ^ as: 


This complex 




i 






(5J 


Thlb form of approximating Q was first suggested by R.T. Jones [1] for 
the case of two-dimensional flow. The matrix Aq can be looked upon as 
representing aerodjmamic stiffness while A- represents aerodynamic damping 
and At represents added mass due to the aerodynamics. The term « 

is an approximation for the time delays inherent in the unsteady aero- 
dynamics and the values of are chosen from the range of reduced fre- 
quencies for which Q(k) has been computed so as to minimize the error in 
this approximation. 

Define the error matrix E(/) as; 


Then the matrices Aq, A 2 ’, “ l,.../are computed so as to 

minimize the spectral norm of ( |)E(-4)1I 2 )* The spectral norm of a 

matrix equals its maximum singular value [,6|]; i.e., || E(iA)|| 2 “ ^ . 

Therefore the maximum singular value of E(f3) is a measure of the size of 

the error in our approximation of Q(>^) and we must choose the values of 

^ such that the corresponding ?'rE(-i)l is the smallest possible, 
m 

The matrices Aq, Aj^, A 2 , and D^, m = 1,.../ are real and are com- 
puted by a least square fit of the aerodynamic data. This is carried 
out as follows [3]; 

Consider the (^,q) element of the matrices, then for this elanent 
we can write (2) as; 


3 


ORIGJNAI PA©E IS 
OF POOR QUALITY 


[I 


cs 

2V 


IMT 




■— — 1 


Substituting for J ■ jk, where J Sot: 


i±Lc 1 


"'i./v 

fit 

iiif 
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Writing out equations like (4) for every k^, i 
bining in a matrix form we get: 

ii^t^ 


l,...n and then com- 
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Since we want the least square solution to (5) to be real valued we 
can write (5) as: 
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The least squares solution to (6) can be obtained using any one of 
the many standard techniques. The program used for the purposes of this 
report Is SNVDEC, developed by NASA^S^ which uses singular value decom- 
position to get Che least squares solution. An explanation of this pro- 
cedure for solving the least squares problem can be found In (6^ . Once 


we have computed A^, A^ 


as the solution to (6), sub- 


‘0’ "1* "2’ m 

stltutlng for Q(/© from (2) In (1) we get the equations of motion to be; 

t 5.-< 




< j L? c, j 


ly) 


From this we can now obtain a state-space model using any of the mini- 
mal realization techniques as 




Ax ^3ii-^rV. 



5 


whflfiC X " atatt vactor 

u ■ control vactor 
2 ■ guat vactor 
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MODIFICATION FOR RIGID BODY ANALYSIS: 

Slnca tha fraquanciaa for tha rigid body motion of tha vahlcla ara 
small comparad to chat of tha flaxura modas, It la Important Co gat a good 
fit for Q(^) at low fraquanciaa In ordar to analyza tha rigid body motion. 

Ona modification In tha abova procadura, which would achlava this, la 
Co aat Aq " Qf((0) so that wa now naad Co obtain a bast laaat squaraa fit 
for only Aj^, A^, D^, m - 1, g . 

Using this approximation for A^, aquation (6) la modlfiad aa: 
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r im^c. 


Ay,,v 
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Thm fit for Aj^, A^, D^, m " 1,...^ Chun obtained will give a good 
approximation for tha rigid body dynaolca of the aircraft. 


MODEL FOR PAST ARW-2: 

The DAST (Dronea for Aerodynamic and Structural Teatlng) ARW-2 
(Aeroelaatlc Reaearch Wing - Number 2) haa a high aapect ratio (10.3) 
atipercrltlcal wing with a 25 degree sweep at midchord mounted on a Fire- 
bee drone fuaelage. 

The unateady aerodynamic Influence coefficient matrix for DAST ARW-2 
at a flight condition of M * 0.86, M “ 15000 ft was provided by NASA - 
Che model consisted of two rigid body modes (plunge and pitch), 10 symmetrl'' 
elastic modes, three control surfaces (stabilizer, outboard aileron. In- 
board aileron) and one gust state. 

The aerodynamic forces were computed using an aerodynamic/ structural 
Interface and a doublet lattice aerodynamics code, contained In the ISAC 
program, for twelve reduced frequencies (0.0, 0.05, 0.1, 0.2, 0.3, O.A, 

0.5, 0.6, 0.7, 0.8, 1.0, and 1.2). 

The elemenxs of Che aerodynamic coefficient matrix were plotted on 
polar plots with the magnitude In decibels versus the p ^e. Several of 
these plots which are typical of Che first column of Q(^) (the force due 
Co plunge) are attached in Fig. 1. Plots typical of forces due Co the 
other modfcs and Che control surface deflections are shown In Fig, 2. 

From Che plots In Fig. 1 It was seen Chat the forces due Co plunge 
are very small at low frequencies compared Co Che forces at higher fre- 
quencies. Therefore the fit for the first column of Q(<i) i4s obtained 
using Che rigid body modifications described In Che earlier section, l.e., 
the first column of Aq (see Eq. (2)) was forced co be equal to the first 
column of Q(-^ at k • 0. 
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The fits for QCi) ware obtained for different values of and a comparison 
of the maximum singular values of the difference between the actual values 

A A 

of Q(d) and the approximations for different k^ showed that using k " 0.13 
gives the best approximation. The following table lists the values of the 

A 

maximum singular value at zero reduced frequency for some values of k. 


S.N. 

k 

a[E(0)I 

1 

0.1 

214.4 

2 

0.125 

138.2 

3 

0.13 

137.0 

4 

0.135 

147.0 

5 

0.15 

162.8 


The approximations of the aerodynamic forces using this value of k 
for the elements plotted In Figs. 1 and 2 are shown In Figs. 3 and 4 re- 
spectively. These plots show that the approximated values follow the actual 
values very closely. Also from these plots we notice that It Is not worth- 
while to use more than one value of k as that would not Improve the fit 
much and would only lead to an Increase In the number of states In the 
state space realization of the model. 

A listing of the program written to obtain the desired fit Is attached 
In the appendix. The matrices A^, A^, A 2 « and corresponding to this fit 
are also listed. Finally an explanation of the parameters of the subroutine 
SNVDEC which was used to obtain these fits Is listed. 


RIGID BODY ANALYSIS FOR DAST AKW-2 


At low frequencies the forces due to the flexure of the aircraft are 
small and the equations for the rigid body motion (with the controls fixed) 
can be written as: 
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T 
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where M, A^, A^, A 2 , are 2x2 matrices corresponding to the rigid body 
modes and (plunge and pitch) . 


Defining lag state ^ as ^(s) 




J 



x(-») 


we have 



g 4-x 


Let 
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Then a state space representation of (8) Is given as: 



For the DAST ARW-2 at a flight condition of M ■ 0.86, H ■ 15000 ft we 

have 

c ■ 23.47 In 
q ■ 4.29 psi 
V ■ 908.79 ft/sec 


From the data provided by NASA and the fits obtained, the matrices 
M, Aq, Aj^,‘ A 2 , and Lj for ■ 0.13 are as follows: 




A.-' 


r o*| 07 s’ 7 Xl<^' 0-(30^n«l0 

o-n73t<io* 


Ax* 


0-376C3SfXlO* o.|3a«SXKJ^ 
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0780027X^0 

-O-lSIf/*?’*'®' 






0 ).' 


0dii5'<07 yiO 
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25 


Using these values, a state space realization as In (9) was obtained 
for the DAST ARW-2. 

The roots of the system for this realization were found to be 
6.1146, 1.828 X 10“^, -5.1767 + J7.5443, -1.1977 x 10^ + JO. 17228. 

The last two roots correspond to the lag terms and the other complex 
pair to the short period motion of the aircraft. 

Neglecting the lag terms, l.e., taking only the top left 4x4 part 

_2 

of the matrix in (9), the roots of the system were found to be 3.391 x 10 , 

7.5875 X 10~^, -1.1188 + J3.5572. 

A possible explanation for the discrepancy between the roots for the 
rigid body motion with and without the lag terms Is that the lag terms be- 
come Important only at high frequencies and at those frequencies the flex- 
ure of the aircraft Is large and we cannot neglect the forces due to the 
flexure modes. 

CO NCLUSION 

A procedure for modelling unsteady aerodynamics of a flexible vehicle 
for t O' ■ '"ol system design was developed. Using this procedure a model for 
DAST-ARW-2 was obtained and this model was found to approximate the actual 
unsteady aerodynamic forces very closely. 

Using this model the eigenvalues for the rigid body motion of the 
aircraft were calculated and these were found to have good correspondence 
with those calculated by NASA. 

This model can now be used for designing active flutter suppression 
control systems for the DAST ARW-2. 
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-.22879SI8E«01 
.9794I292£*0I 
-.84201844E«0I 
-.183S2327E«02 
.I0427999C602 
.23289894C«02 
.24724347E602 
-.97438098E602 


-.23712S34£»03 
-.I4588059E603 
.43269445E«01 
-. 108998 1 4E603 
.I94I4482E602 
1 727021 7E+03 
.99993304E601 
-.I0077923E+04 
-.738t1418E*02 
-.:2413784£*03 
-.113318802603 
-ir52l4y44EV04 


.17539497E603 

.82194747E602 

-.48302880E600 

.93393017E601 

.24792928E-01 

.33002:07E600 

-.28334319E602 

-.12418812E602 

-.93898S44E602 

-.I9247093£t03 

.41444428E602 

.424304132601 


. 7793*4212602 

.249397492602 

.333479472601 

.5524777!c*01 

-.43302780E600 

-.734084132601 

-.13099147E602 

.40972982E602 

.I0004S90E6O2 

.47494049E601 

.204010902601 

.477397822602 


.1314S939E604 
-.9>840S4.3E602 
.29319449E603 
-.21919444E603 
-.188937342602 
-.28941 1UE 603 
.13894489E603 
.11777249E604 
-.313089992602 
.71440542E601 
.94440807g6Q? 
'.l8o94“l05E604 


.34578743E603 
.178I3149E603 
.30944404E60I 
. 497479 18E602 
-.3271 401 4E600 
-.27293774E601 
-.10442332E603 
.913229702602 
.487291422602 
-.92787777E600 
.473912872601 
.11 29837 2E602 
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ORIC'.NAl PAGE J3 
OF POOR QUALITY 

J ■ . • 

' • » , 

-.505l9gA0€*03 .70I40330?»0' -.»029I392S+0J .:324;^*2t+03 

- . 1 4«73790e^03 . < 04l4434c*02 1 4478<>47c*0i . 4doJ4'. 74c-r02 

-.533274I4£*03 .821 758t 9£402 •.1038840«E*03 .9907«)90. P.402 

-.I3293320E403 .28707300£402 -,42219490€»02 .«?31 : 212E40I 

k337449i8E402 •.30897333c40t . I 2340S«7£40l .U37i328€40t 

-.U742933E«03 .32034933E402 •.2535U67E402 -.79749038I-402 

-.U40823w£402 . t8780049E402 U734H2E402 -.t;410904E402 

.779A7S07Et03 .3A2S8390E«03 -.2SS43A2IE403 .47l47343E-»03 

•.7S729487E40t -.40492408E44t .42284S4IE40I .479983 19E 40 I 

.3998703tE402 .t4470247«402 .1801 1440E402 •.I444d380<i402 

.980223S9E402 -.AtlUSOOEiOl .43931740E402 -.22H3728E402 
-.23I73931E403 -.3220401 lEtOl .32097028E403 -.1017IC23E403 

i* 

<«•«•••••••• 15TM C0co4h •«•«•••««•«•••*• 

■t 

-.342434I4E404 .4570794IE*Oi 77473 lAOE+OI .22398733E404 

.42874289E442 .I4348737E403 -.425113448403 -,2877t909£403 

-.447444828403 -.818923848402 .449283148402 .4740«712e403 

.713732998403 .482332478402 -.834427208402 -.317304028-403 

.438708328402 -.327377378400 .241 231448401 -.308240338402 

.888144508403 -.443417758401 .944303418402 -.437214078403 

.808444178402 -.133723578403 .325798478403 .199235058403 

.107030248404 -.833434008403 .123397098404 . l04o7;08E404 

.120228818403 .444833958402 -. i 004 18838403 -.199471448-403 

-.434740278403 .410748918402 -.144178118403 .442590138402 

-.A35543448403 .338252948402 -.943188948402 .204430448-403 

.216219528404 -.944951088403 .141341018404 .223397108404 

:« 

:M«6666«46 14TM COLUNM 

.222433488403 .297879038404 -.320949348403 -.244144328405 

.445818808403 .805734228403 -.744448108403 -.114050298404 

.499997778404 .223438418404 -.177539948404 -.723825248404 

-.181182208404 .742851838403 -.134118398404 .I291960«8404 

-.243383088403 -.247703338402 -.933174748400 .31 T7lS45€403 

-.314190898404 -.393172128403 -.989727388403 .393842338404 

.100072528403 -.140937998403 -.397854498402 .414228148402 

.128328048403 .747414048404 -.830883898404 -.200849458405 

.142902378403 .501 144718402 -.491800478402 -.218095378403 

.433345388403 -.529789498403 .483037038403 -.437475978401 

.111209708404 -.447841058403 . 748241 658403 -.704418548403 

.158033028403 -.118235278404 .53. 724538404 -.14.523498405 
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Dtaoriptlon ; Th« purpoM of SNVDEC is to oonputs ths slngulsr-vslus dsooaposl* 
tion (rs^. 2) of a raal ■ x n (a 2 n) aatrix A by parforaing tha faotorl- 
lation, 


A ■ UQV* 


whara U is an a x n aatrix whosa ooluana ara n orthonoraalizad aigan- 
vaotora aasooiatad with tha n largast aiganvaluas of AA', V is an 
n X n aatrix whoso ooluana ara tha orthonoraalizad aiganvaotora assooiatad 
with tha n aiganvaluas of A' A, and 


Q ■ dlag 


whara a* (1 ■ 1,2,... ,n) ara tha nonnagatlva squara roots of tha aigan- 
valuas or A* A, oallad tha singular valuas of A. Options ara provldad for 
tha ooaputation of rank A, slnipilar valuas of A, «i orthonorsal basis for 
tha null spaoa of A, tha psoudoinvarsa of A, and tha laast squaras solution 
to 


a ■ B 


Both A and B ara storad as varisbla-dlBansionad two-dlsansional arrays. 

Tha oooputational prooadura is dasoribad in rafaranoa 2 on psgas 135-151. 
Basioaliy, Housaholdar transforaations ara appllad to raduoa A to bidlagonal 
fora aftar whioh a OR algoritha is usad to find tha singular valuas of tha 
rsduoad aatrix. Combining ra suits givas tha raquirad oonstruotion . 

Souroa of ^ftwa^ t LaRC Analysis and Computation Division subprogrsa library 
with Bodificatlons by Bmast S. Arnatrong, LaRC 

Calling saquanoa : CALL SNVDBC(IOP,MD,IID,M,N,A,NOS,B,IAC,ZTEST,Q,V,IRANK, 

APLUS,lEm) 

Input arguaants t 

ZOP Option ooda: 

1 Tha rank and singular valuas of A will bs ratumad. 

2 Tha aatrioas U and V will ba ratumad in addition to tha 

inforsation for lOP ■ 1. 

3 Xn addition to tha inforsation for lOP ■ 2, tha laast 

squaras solution to AX ■ B will ba ratumad. 

4 Tha psaudoinvarsa of A will ba ratumad in addition to tha 
inforsation for lOP ■ 2. 

Tha laast squares solution will ba ratumad in addition to 
tha inforsation for lOP « 4. 
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MD The maximum firat dimension of the array A as given in the 

DIMENSION statement of the calling program 

ND Maximum first dimension of the array V 

M The number of rows of A 

N The number of columns of A 

A Matrix stored as a variable-dimensioned two-dimensional array. 

Input A Is destroyed. 

NUS The number of column vectors of the matrix B 

B Two-dimensional array that must have row dimension at least NOS 

in the calling program. B contains the right aides of the 
equation to be solved for lOP s 3 or lOP s 5. B need not 
be input for other options but must appear in the calling 
sequence . 

lAC The number of decimal digits of accuracy in the elements of the 

matrix A. This parameter is used in the test to determine zero 
singular values and thereby the rank of A. 

Output arguments : 

A On normal return, A contains the orthogonal matrix U except 

when lOP s 1. 

B On normal return, B contains the least squares solution for 

lOP a 3 or lOP a 5. 

ZTEST The zero test computed as ||A|| * using the matrix 

Euclidean norm except when N a 1 . When N a 1 , 


ZTEST a 


Q A one-dimensional array of dimension at least N which upon return 

contains the singular values in descending order 

V A two-dimensional array that must have first dimension ND and 

suoond dimension at least N. Upon normal return, this array 
contains the orthogonal matrix V except when lOP s 1. The 
last N - IRANK columns of V form a basis for the null space 
of A. 

IRANK Rank of the matrix A determined as the number of nonzero singular 
values using ZTEST 


ORlG!f,VAL PAGI fs 
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APLUS 


lERR 


A two-dimensionai array or first dimension ND and second dimension 
at least M. Upon normal return, this array contains the pseudo- 
inverse of the matrix A. If lOP does not equal 4 or 5, this 
array need not be dimensioned, but a dummy parameter must appear 
in the calling sequence. 

Error indicator; 

TERR =0 A normal return 

lERR s K •' 0 The Kth singular value '.las not been found after 
30 iterations of the QR algorithm procedure. 
lEKR = -1 Using the given lAC, A is close to a matrix 

which is of lower rank and if the aoouraoy is 
re<*uoed, the rank of the matrix may also be 
reduced . 


CWMOM blocks ! None 

Error messages ; None. The user should examine lERR after return. 

Field length ; 2072 octal words (1082 decimal ) 

Subroutines employed by SNVDEC ; None 

Subrou tires employing SNVDEC : FACTOR, CTROL, CSTAB, Dc^'AB, DISCREG 

Conaaents ; SNVDEC may be applied to matrices stored as one-dlmenslonal arrays 
by setting MD s M and ND s N In the calling sequence. 

The subroutine Is Internally restrlct,ed to N $ 150. 
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Figure 1.- Vertical coeponent of clastic mode shapes. 


38 


O!! F/'Cl.- 
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(b) «»)« 2. frequency - U.7 he. jeneraliied »ss • 0.526 Ib-sec /in 



•ode 3. 




Figure 1.- Ccntlnued. 
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(f) mode 6, frequency * 44. E hr, generalized mass » 0.593 lb-sec‘/in 
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Figure 1.- Continued. 
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ORIGINAL FU iJ. 1 
OF POOR QUALliY 



Ch) ■ode 8. frequency - 66.6 hr, generalized wiss =• 0.0621 Ib-»ec /in 


^#1 
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(j) node 10, frequency - 80.3 hi. gencraliicd »ass - 0.696 lb-sec‘/in 
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Wo .’.king Pap«r No. 5 
Modified Outboard Control 
Surface Transfer Function 


W. I,. Garrard 
Department of Aerospace 
Engineering emd Mechanics 
University of Minnesota 
Minneapolis, Minnesota 


This paper discusses the modified transfer function for the 
outboard aileron. The new transfer function is 

xio / „ xio'^( f za.6s ) 

Note that the numerator term and two of the denominator terms 
have the same nondamped natural frequency, the deunping factor of 
the numerator is 0.03 while the damping factors associated with 
the denominator terms of the same frequency are 0.23 and 0.3. 

'^hus the numerator dynamics effectively cemcels the phase shift 
resulting from the one of the denominator factorn and except at 
near the natural frequency of 477.5 rad/sec the numerator dynamics 
also cancel the 40 db/decade roll-off resulting from one of 

the denominator factors. Bode plots, for the response of the exact 
model of the outboard aileron are shown in Figs 1 and 2. 



OF POOR QUALITY 

The most obvious approximation is the cancellation the 
numerator dynamics with the most lightly d^unped of the quadratic 
terms of the same frequency in the aenominator. This results in 
the 5th ordar model below 

- 4 , 0 -^ )(/0 

C^-^J8oXa^ ^ZS\S 28(iS4(^n,5f) 

The Bode plots for this transfer function are gi’ en in Figs 3 
\nd 4. Phase is almost the same as for the exact model and gain 
only differs near 477.5 rad/s. A. 3rd order appro iximation results 
from neglecting the highest frequency term in the denominator. 

This results in the approximation 

T 77 '^ /O ^ 

(s-^iSo)Cs^ -hzsls 

Bode plots for this transfer function are shown in Figs 5 and 6. 

Up to approximately 300 rad/ sec this gain 6md phase are the 
same as for the exact model. Since the 6th structural mode has 
a frequency of 225 rad/s at the flutter frequency it is felt that 
the 3rd order actuator model is adequate for initial control 
studies. 
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